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ABSTRACT 

This paper represents the investigation on the response of symmetric laminated composite plate. In this paper, 
an attempt is made to study the Buckling characteristics of a laminated composite plates, using layer wise HSDT. 
Buckling is one of the major modes of failure of laminated structures and it is necessary to predict the critical load of the 
structural component, for the easy replacement with good load bearing capacity. Studying buckling characteristics of 
laminated composite plates is essential, in order to have a deep insight of their behaviour for proper analysis and design. 
The derivation of equations of motion for layer wise higher order models is obtained, using the virtual work principle. 
Further, the analytical solutions are developed to find buckling characteristics of laminated composite plates. The 
solution predicts critical buckling load of laminated composite plate. The results obtained by using the present model are 
compared with the available literature. 
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1. INTRODUCTION 

Laminated composites with continuous fibres are widely being used in various engineering fields, such as 
aeronautical and aerospace industry, marine, aviation, civil, sport, as well as in other fields of modern technology 
and other applications. They are preferred due to their characteristics like, high stiffness to weight ratio, excellent 
fatigue strength, high energy absorption, self damping capacity, high strength to weight ratio, low weight, 
resistances to electrochemical corrosion and good resistance to corrosive agents, capable of being engineered 
according to requirements. It is a challenging task, to find the accurate prediction of the response characteristics of 
composite structures. Hence, it is necessary to analyze the buckling characteristics of laminated composite plates. 
In the present paper, the displacement model is developed using the Layer wise HSDT. The equations of motion 
are derived using the principle of virtual work. The solutions are obtained using Navier’s method. Later buckling 
characteristics are determined using Layer wise HSDT. Linally, the accuracy of the proposed theory is verified 
based on the results available. 

2. THEORETICAL FORMULATION AND DEVELOPMENT OF LAYERWISE HSDT 

There has been a considerable progress in understanding behaviour of composite laminates. It is noticed 
that, anisotropic multilayered structures posses transverse discontinuous mechanical properties and higher 
transverse shear and transverse normal stress deformability. In order to model materials with such behaviour, two 
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different approaches have aroused, they are equivalent single-layer theories (ESL) and layerwise theories (LWT). In order 
to overcome the limitations of CLPT and FSDT, Higher-order Shear Deformation Theories (HSDT), which involve higher- 
order terms in Taylor's expansion of the displacements in the thickness coordinate were developed. These models can 
disregard shear correction factors and give more accurate and stable transverse shear stresses. It is noticed that, LWT 
models have some advantages over the conventional 3D models. First, they allow independent in-plane and through the 
thickness interpolation, the element stiffness matrix can be computed much faster. Second, even the volume of input data is 
reduced, LWT are capable of achieving the same level of solution accuracy, as a conventional 3D models. 



Figure 2.1: Layer Wise Kinematics of ‘N’ Layer Laminated Composite Plate 


A rectangular plate having width a, length b and thickness h along x, y and z axes, respectively is considered. 
Figure 2.1, shows the one dimensional layer wise kinematics of three layered composite. The displacement 
components t/ 2 \ v^ and w^ 2 \ along x, y and z directions, respectively for the middle layer are expanded using, Taylor’s 
series in terms of thickness coordinate z^ 2 \ 


Displacement Field 

u (k) (x,y,z) = u 0 {x,y) + z w 9 x (x,y) + (z (k) ) 2 u 0 *O,y) + (z^f 9 x *(x,y) + 0 X 

u®(x,y,z) = v 0 (x, y) + z^ k) 9 y (x,y) + (z (k) ) 2 iVfey) + (z (k) ) 3 0 y *(t:,3') + 
wW(i:,y,z) = w 0 (x,y') 

( (Sr=middie layer T 0 [m)) for layers above middle layer 
Where ® m = \ .... , 2 .. 

(S'i mi e ayer ~ 0 (m)) for layers below middle layer 


Where u 0 , v 0 and w 0 are the displacements of the middle plane along x, y and z directions, respectively. The 
parameters 0 X and 9 y are the rotations of normal to middle plane of k lh layer, about x and y axes, respectively (k=l, 2, 3...,) 
and u 0 *, v 0 *, 8 X * and Q y are the higher-order terms in the Taylor’s series expansion and represent higher-order transverse 
cross sectional deformation modes, of the middle layer. Laminated plate is made of laminae, having a fibres oriented at an 
angle 0, measured from the material x to global x axis. The stress strain relations of the laminae are therefore, defined in 
material coordinate system as: 
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In which 

o= ((O!, o 2 , t 12 , t 13 , ^ 23 )’ are the stress components of k th laminae in material coordinates 
e = ((£j, £,, y 12 , Y 13 , 723 )*) are the strain vectors k " 1 laminae in material coordinates and 
Qj’s are matrix of material elastic coefficients for k lh laminae, given as, 

C — ^1 C — ^ 21^1 


c 


22 


E 2 

(l — M-12M-2l) 


C 3 3 = G 12 ;C44=G 2 3;C 5 5 = G 13 

Since, all quantities should be referred to a single coordinate system; it is needed to establish transformation 
relations among stresses and strains, in global system to the corresponding quantities in material (local) coordinate system. 
The constitutive matrix in global coordinate system will then be of the form: 
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Where, 

Qii=C„*Cos[0] 4 +2*C 12 *Sin[0] 2 *Cos[0] 2 +4*C 33 *Sin[0] 2 *Cos[0] 2 +C 22 *Sin[0] 4 ; 
Q 12 = (Cn+C 22 -4*C 33 )*Sin[0] 2 *Cos[0] +C 12 (Sin[0] 4 +Cos[0] 4 ); 

Q 22 =C„*Sin[0] 4 + (2*C 12 +4*C 33 )*Sin[0] 2 *Cos[0] 2 +C 22 *Cos[0] 4 ; 

Qu= (C u -C 12 -2*C 66 )*Sin[0]*Cos[0] 3 +(C 12 -C 22 +2*C 33 )Sin[0] 3 *Cos[0]; 

Q 23 = (C 11 -C 12 -2*C 33 )*Sin[0] 3 *Cos[0]+(C 12 -C 22 +2*C 33 )Sin[0]*Cos[0] 3 ; 

Q 33 = (Ci 1 +C 22 -2*Ci 2 -2*C 33 )Sin[0] 2 *Cos[0] 2 +C 33 (Sin[0] 4 +Cos[0] 4 ); 
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Q 4 4=C44*Cos[0] 2 +C 5 5*Sin[0] 2 ; 

Q45=(C 55 -G w )*Cos[0]*Sin[0]; 

Q 5 5=C55*Cos[0] 2 +C 4 4 :|< Sin[0] 2 ; 

3. BUCKLING ANALYSIS OF LAMINATED COMPOSITE PLATES BASED ON LAYERWISE 
HSDT 

Buckling in structures is due to in-plane loads, represented by mechanical and thermal loads. Buckling is one of 
the major modes of failure of laminated structures and it is necessary to predict the critical load of the structural 
component, for the easy replacement with good load bearing capacity. The objective of buckling analysis, is to determine 
the critical buckling load. The buckled configuration exhibits the structure, when the in-plane load exceeds its critical 
buckling load and leads to instability. Hence, the effect of thickness ratio, aspect ratio and modulus ratio on the buckling of 
laminated composite plates are investigated 

4. NAVIER SOLUTION USING HIGHER-ORDER DISPLACEMENT MODEL BASED ON 
LAYERWISE THEORY 

In the Navier method, the displacements are expanded in a double Fourier series, in terms of unknown parameters. 
The choice of the trigonometric functions in the series is restricted to those which satisfy the boundary conditions of the 
problem. Substitution of the displacement expansions in the governing equations, results in an invertible set of algebraic 
equations, among the parameters of the displacement expansion. 

The simply supported boundary conditions for the higher-order shear deformation theory are: 

At edges x = 0 and x = a 

v 0 = 0, w 0 = 0, 0 y = 0, M x = 0, v 0 = 0, 0 y * = 0, M x * = 0, N x = 0, N X T = 0, 

At edges y = 0 and y = b 

u 0 = 0, w G = 0, 0 X = 0, M y = 0, u 0 * = 0, 0 X * = 0, M y = 0, N y = 0, N y * = 0, 

The simply supported boundary conditions shown in equations are considered for solutions of laminated 
composite plates, using displacement model. The boundary conditions are satisfied by the expansions: 

oo oo 

u Q ( x,y,t) = II U mn (t) cos ax sin fiy 

m=l n= 1 

v 0 (x, y, t) = II V mn (t) sin Of cos/?v 

m= 1 n= 1 

w 0 (x,y,t) =11 W mn (t) sinm- sin/?v 

m=l /(=! 
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d x (x, y,t)= EE X mn (t) cosax sin py 

m—\ n =1 

B y ( x,y,t) = EE Y mn (t) sinm: cos py 

m =1 n =1 

Uq ( x,y,t) = EE U mn (t) cos ax sin py 

m =1 n =1 

V a *(x,y,t) = EE Vl(0 sinm- cos/?y 

m=l n= 1 

B* (x,y,t) = EE Y*„(0 cos a* sin/A 

m=l n=l 

The mechanical loads are also expanded in double Fourier sine series as: 


q (. X,y,t ) = EE £L,(0 sin ax sin py 

m= 1 n=l 


Where Q mn (z, t) = 


4 

ab 
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j j q(x,y,t) 


0 0 


sin (Xx sin jBy dxdy Where a = 


ni7T 


a 


and (3 = 


n 71 

b 


The governing equations of motion are obtained as: 

(A n a 2 + A^ 2 )u mn +(A 12 + A 33 )V mn a/3 + (B u a 2 +B 33 /3 2 )x mn +{B l2 +B 33 )Y mn a/3 + 
(A H a 2 + A 36 j3 2 )u:„ + (A IS + A 36 )v: mn aj3 + (B u a 2 + B 36 p 2 )x* mn + (b 15 + B 36 )C aj3 = 
(hu^ + hx^ + i.iil + i.rJ-iN 1 ) 

(A 2 1 + A 33 ) U mn a(3 + (a 33 oc~ + A 22 /A ) V mn + (B 21 + S 33 ) X mn ap + (b 33 or + B 22 p~)Y nn + 

(A 24 + Afi )U + (^36°^ +A 250 2 ) V m n +( B 24 + S 36 a P + ( 5 36 + S 2s/A I 1 '™, = 

fc, +^t n +/ 3 C+l4Yl)-(iV 2 ) ) 

tn « 2 +L 22 ^)w mn+ L 11 «X mn +L 22 pY mn +2L 13 «f/l + 

2L 24 J3V:+3L I5 <„,+3L 26 pY* mn = {pW mn +<2 m „ ) 

(BjjQT 2 + B 33 P 2 )u mn +(B 2 i + A 33 )V m n ^ + L 11 W mn a+(D ll « 2 +D 33 ^ 2 +L U )x mn + 
(Dn + £> 33 )A„ cc/5 + (s 41 or 2 + A ,/? 2 + 2 L, 3 + {B 5l + 5 63 )<„ a/3 + 

(D u a 2 + D }6 p 2 + 3 L 15 )z, 1 „ + (d 15 + Djr a(3 = 

M™+a +/ 4 a;„+ 1 5 x:„)- (m 1 ) 


www.tjprc.org 


editor@tjprc.org 





68 


S. Lokesh, GS Vivek & L. Praveen 


{B n + B 33 )U nm afi + (b 33 a- + B 22 fl- )V mn + Y 22 W mn /3 + (D 2 | +D 33 )X mn a/3+ 

(d 33 a- + d 21 p- + L 22 )y mn + (s 42 + B 63 )u mn afi + 

K<* 2 +B 52 j3 2 +2L 24 )v; m +{D 36 +Djxlap + 

(D 36 a 2 + IXjf- + 3 L 26 )C = (r 2 V mn +1 3 Y mn +1 4 V; nn + 1 5 Yl)- (M 2 ) 

(A 4l a 2 +A 63 p 2 )u mn +{A 42 +A 63 )v mn ap + 2L 3l W mn a + (b 41 « 2 + 5 63 /? 2 + 2L 3l )x mn + 

( B 42 + B 63 ) Y mn a P + (Am#" + A 66^~ + 4 ^33 )^ mn + ( A 45 + A 66 ) V mn a P + 

(b 44 cc 2 +B 66 p 2 + 6 L 35 )x* nm + (fi 45 + fl 66 )C = 

M™+j 4 x nm+ i 5 u: n+ i 6 x: n )- { N-) 

Ust + A 63 ) U mn afi + (a 63 « 2 + A 52 fi 2 ) V mn + 2 L 42 W m J + (B S1 + B 63 ) Y m „ or/? + 

(b 63 +B 52 fi~ + 2L 42 )Y mn + (A 54 + A 66 )u mn aft + 

(A 66 er + A 55 /?- + 4L 44 )v mn + (S 54 + B 66 )X mn af3 + 

(B 66 a 2 + B 55 p 2 + 6L 46 )C = (l 3 V mn +1 4 Y mn + / 5 Vi +1 6 Yl) - (N* 2 ) 

[b u u 2 + B 36 fi 2 ) U mn + (B 24 + 5 36 )V mn afi + 3L 5l W mn a + (D 4l a 2 +D 63 p 2 + 3L 5l )x mn + 

(D 42 + D 63 )Y mn aft + (B u a 2 + B,J 2 + 6L 53 ) u; im + (S 54 + S 66 )vj„ afi + 

K« 2 +d 66 /? 2 +9l 55 )x; nn +(D 45 +D 66 )C a(3 = 

(*15 + * 36 ^™ afi + ( B 36 «r 2 + B 25 ^ 2 )v mn +3L 62 W,„,^+(D 51 + D 63 )X m „ 00 + 

(d 63 a 2 + D 52 ff- + 3L 62 )f„ + (fi 45 + S 66 )£/>/? + 

(B bb a 2 + S 55 /? 2 + 6L 64 )v; n + (D 54 + D 66 )x* mn a/3 + 

K a 2 + D 55 p 2 + 9 l 66 )C = (/ 4 y,„„ + / 5 Y mn + / 6 y,: +i 7 y;„ ) - (m’ 2 ) 

The analytical solutions for buckling problems of laminated composite plates, can be obtained by setting the 
distributed load over the surface of the plate and time derivatives to zero: 
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Where, K = N cr {2 x a 2 +i 2 /? 2 ) 

The in-plane loading conditions for laminated composite plates are given as: 

Uni-axial compression (UAC): 2 X = — 1,2 2 =0 
Biaxial compression (BAC): 2 X = —1,2-, = — 1 

The elements Sijand (i=l,2....9 and j=l,2_9) are given and solutions for the equations for each m,n gives U mn , 

v m „, w mn , X mn , Y mn , u mn , V inn , X mn , Y mn Which are used to compute U o ,V o ,W o ,0 x , 0 y ,U o ,V Q ,0 X ,0 y 

Sn = An*a 2 + A 66 *|3-; S 12 = (Ai 2 +A 6 s)*a *P; 

S 13 = 0; S 14 = B„*a 2 + B 66 *p 2 ; S ls = (B 12 +B 6S )*a *P; 

Si6= A 14 *a' + A36*P“; Sn= (Ai 3 +A36)*a *P; 

S lg = B 14 *a 2 + B 36 *p 2 ; S 19 = (B 15 +B 36 )*a *P; 

S 2 i= (A 2 i+A 66 )*a *P; S 22 = A 66 *cr + A 22 *p~; 

S 23 = 0; S 24 = (B 12 +B 66 )* a *P; 

S 25 = B 66 *a 2 + B 22 *p 2 ; S 26 = (A 24 +A 36 )*a *P; 

S 2 7= A 36 *a' + A 2 5 *P“; S 2 g= (B 24 +B 3 g)^a *P; 
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S 29 = B 36 a 2 + B 25 *p 2 ; S 31 = 0; S 32 = 0; 

S 33 = L,i*a 2 + L2 2 *p 2 ; S 34 = L n *a; 

S 3 5= L 22 *(3; S 36 = 2*L 13 *a; 

S 37 = 2*L 24 * (3; 

S 38 = 3*L 15 *a; S 39 = 3*1^6* P; 

S 41 = B„% 2 + B 66 *p 2 ; S 42 = (B 21 +B 66 )*a *P; 

S 43 = L„*a; S 44 = D n *a 2 + D 33 *p 2 + L„; 

S 45 = (D 12 +D 33 )*a *P; S 46 = B 41 *a 2 + B 63 *p 2 + 2L 13 ; 
S 47 = (B 51 +B 63 )*a *P; S 48 = D 14 *a 2 + D 36 *p 2 + 3L 15 ; 
S 49 = (D 15 +D 36 )*a *P; 

S 5 i= (B 12 +B 66 )*a *P; S 52 = B 66 *a 2 + B 22 *p 2 ; 

S 53 = L2 2 *P; S 54 = (D 21 +B 33 )*a *P; 

S 55 = D 66 *a 2 + D 22 *p 2 + L 22 ; S 56 = (B 42 +B 63 )*a *P; 
S 57 = B„*a 2 + B 52 *p 2 + 2L 24 ; S 58 = (D 24 +B 36 )*a *P; 
S5 9 = D 36 *a 2 + D 25 *p 2 + 3L2 6 ; 

Sgi = A 4 i*a 2 + A 63 *P~; S6 2 = (A 42 +A 63 )*a *P; 

S 63 = 2*L 31 *a; S 64 = B 41 *a 2 + B 63 *p 2 + 2*L 31 ; 

S 65 = (B 42 +B 63 )*a *P; S 66 = A, 4 a 2 + A 66 *p 2 + 4*L 33 ; 
S 67 = (A 45 +A 66 )*a *P; S 68 = B 44 a 2 + B 66 *p 2 + 6*L 35 ; 
S 69 = (B 45 +B 66 )*a *P; 

S 7 i = (A 5 i+A 63 )*a *P; S 72 = A 63 *cr + A 52 *p~; 

S 73 = 2*L 42 * p ; S 74 = (B 51 +B 63 )*a *P; 

S 75 = B 63 *a 2 + B 52 *P 2 + 2L 42 ; S 7s = (A 54 +A 66 )* a *P; 
S 77 = A 66 *a 2 + A 55 *p 2 + 4*L 44 ; S 78 = (B 54 +B 66 )*a *P; 
S 79 = B 66 *a 2 + B 55 *p 2 + 6*L 46 ; 

Sgi = B 14 *a 2 + B 36 *p 2 ; S 82 = (B 24 +B 36 )*a *P; 

S 83 = 3*L 15 *a; S 84 = D 41 *a 2 + D 63 *p 2 + 3L 51 ; 

S 85 = (D 42 +D 63 )*a*P; S 86 = B 44 a 2 + B 66 *p 2 + 6*L 53 ; 
S 87 = (B 54 +B 66 )*a *P; S 88 = D 44 a 2 + D 66 *p 2 + 9*L 55 ; 
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S 89 =(D 45 +D 66 )*a *|3; 

S 9 i= (B 15 +B 36 )*a *P; S 92 = B 36 a 2 + B 25 *p 2 ; 

S 93 = 3*L 62 *P; S 94 = (D 51 +D 63 )*a *P; 

S 95 = D 63 *a 2 + D 52 *p 2 + 3L 62 ; S 9s = (B 45 +B 66 )*a *P; 
S 97 = B 66 *a 2 + B 55 *p 2 + 6 *L 46 ; S 98 = (D 45 +D 66 )*a *P; 
S 99 =D 66 *a 2 + D 55 *p 2 + 9*L 66 ; 

5. RESULTS 


The following material properties are considered for laminated composite plates to perform the analysis. 
Material: Graphite Epoxy 

Young's Modulus: E! =25 Gpa,E 2 =lGpa 
Shear Modulus: G i2 =G 23 =Gi 3 =0.5Gpa, 

Poisson's Ratio: fd I2 = /l 23 = fd 13 = 0.25 

The numerical results obtained from the buckling analysis are shown in graphs 



Figure.5.1: Comparison of Non- Dimensional Uni-Axial Buckling Load of a Simply 
Supported Cross-Ply Plate (0 0 /90°/0 0 ) with Different Modular Ratios 



Figure.5.2: Effect of Aspect Ratio and Modular Ratio on Uni-Axial Buckling Load of a Simply 
Supported Cross Ply (0 0 /90 0 /90°/0 0 ) Laminated Composite Plate 
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Figure.5.3: Effect of B/H Ratio and Modular Ratio on Non- Dimensional Uni-Axial Buckling Load of a 
Simply Supported Cross Ply (0°/90 0 /90°/0 0 ) Laminated Composite Plate 



Figure.5.4: Effect of Aspect Ratio, Modular Ratio and A/B Ratio on Non- Dimensional 
Bi-Axial Buckling Load of a Cross-Ply Laminated Composite Plate 



0 20 40 gyjj 60 80 100 120 


Figure.5.5: Effect of Aspect Ratio, on Non- Dimensional Uni-Axial Buckling 
Load of a Angle-Ply Laminated Composite Plate 
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Figure 5.6: Effect of Modular Ratio, On Non- Dimensional Uni-Axial Buckling 
Load of a Angle-Ply Laminated Composite Plate 



Figure5.7: Variation of the Critical Non- Dimensional Uni-Axial Buckling Load with The in-Plane Ply 
Modulus Ratio E x /Ei in the Kange[ 1,40|,A/H=10 With Different Ply Stacking Sequences 

6. CONCLUSIONS 

A Layer wise higher order shear deformation theory, is successfully developed for the Buckling analysis of 
laminated composite plates. Many problems are solved covering different features of laminated composites, such as ply 
orientations aspect ratio, thickness ratio and loading etc. The analysis has verified that, the present LHSDT is capable to 
predict the buckling response of laminated composite plates, as that of the other theories. It has been shown that, the effects 
of side-to-thickness ratio, aspect ratio, modular ratio play an important on non-dimensional buckling loads of laminated 
composite plates. The information derived from the study can help designers in choosing a particular side-to-thickness 
ratio, aspect ratio, modulus ratio in the material for designing a laminated composite plate. 
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